Control of chaotic transport in Hamiltonian systems 
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It is shown that a relevant control of Hamiltonian chaos is possible through suitable small per- 
turbations whose form can be explicitly computed. In particular, it is possible to control (reduce) 
the chaotic diffusion in the phase space of a Hamiltonian system with 1.5 degrees of freedom which 
models the diffusion of charged test particles in a "turbulent" electric field across the confining mag- 
netic field in controlled thermonuclear fusion devices. Though still far from practical applications, 
this result suggests that some strategy to control turbulent transport in magnetized plasmas, in 
particular tokamaks, is conceivable. 

PACS numbers: 05.45.Gg; 05. 45. Ac; 52.25.Xz 



Transport induced by chaotic motion is now a stan- 
dard framework to analyze the properties of numerous 
systems. Since chaos can be harmful in several contexts, 
during the last decade or so, much attention has been 
paid to the so-called topic of chaos control. Here the 
meaning of control is that one aims at reducing or sup- 
pressing chaos inducing a relevant change in the trans- 
port properties, by means of a small perturbation (either 
open-loop or closed-loop control of dissipative systems 
[ilyl) so that the original structure of the system under 
investigation is substantially kept unaltered. Control of 
chaotic transport properties still remains an open issue 
with considerable applications. 

In the case of dissipative systems, an efficient strategy 
of control works by stabilizing unstable periodic orbits, 
where the dynamics is eventually attracted. Hamiltonian 
description of microscopic dynamics usually involves a 
large number of particles. However, methods based on 
targeting and locking to islands of regular motions in 
a "chaotic sea" are of no practical use in control when 
dealing simultaneously with a large number of unknown 
trajectories. Therefore, the only hope seems to look for 
small perturbations, if any, making the system integrable 
or closer to integrable. In what follows we show that it 
is actually possible to control Hamiltonian chaos by pre- 
serving the Hamiltonian structure. 

Chaotic transport of particles advected by a turbu- 
lent electric field with a strong magnetic field is asso- 
ciated with Hamiltonian dynamical systems under the 
E X B guiding center approximation Although it 
has been shown that the E x B drift motion of the so- 
called guiding center can lead to a diffusive transport in 



a fairly good agreement with the experimental counter- 
parts 0, , it is clear that such an analysis is only a first 
step in the investigation and understanding of turbulent 
plasma transport. The control of transport in magneti- 
cally confined plasmas is of major importance in the long 
way to achieve controlled thermonuclear fusion. Two ma- 
jor mechanisms have been proposed for such a turbulent 
transport, transport governed by the fluctuations of the 
magnetic field and transport governed by fluctuations of 
the electric field. There is presently a large consensus to 
consider, at low plasma pressure, that the latter mecha- 
nism agrees with experimental evidence 0. In the area 
of transport of trace impurities, i.e. that are sufBciently 
diluted so as not to modify the electric fleld pattern, the 
present model should be the exact transport model. Even 
for this very restricted case, control of chaotic transport 
would be very relevant for the thermonuclear fusion pro- 
gram. The possibility of reducing and even suppressing 
chaos combined with the empirically found states of im- 
proved confinement in tokamaks, suggest to investigate 
the possibility to devise a strategy of control of chaotic 
transport through some smart perturbation acting at the 
microscopic level of charged particle motions. 

First, we briefly describe the Hamiltonian with 1.5 de- 
grees of freedom modeling the E x B motion of charged 
test particles in a "spatially turbulent" electric fleld. 
Then we formulate the problem of control and analyt- 
ically derive the partial control term for a Hamiltonian 
describing the motion of these test particles. Finally, we 
report the numerical evidence of the effectiveness of the 
method. Let us begin by describing the model whose 
dynamics we want to control. In the guiding center ap- 
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proximation, the equations of motion of charged particles 
in the presence of a strong toroidal magnetic field and of 
a nonstationary electric field are 



dfx 



^ FC^ R - ^ f-dyV{x,y,t) 



(1) 

where V is the electric potential, E — —WV, and B — 
Be^. To define a model we choose 



/ such that H + f is integrable. In this article, we are 
interested in finding a partial control term /2 of order 
such that the Hamiltonian given by He = Hq + eV + 
is closer to intcgrability, i.e. such that He is canonically 
conjugate to Hq + 0{e^). We perform a Lie transform on 
He, generated by a fmiction S: 



H'^ = e'^He = He + e{S, He} + -{S,{S,H,;\} + ...., (6) 



y(x, t) = ^ sin[k ■ X + (^k - c^(k)t]. 



(2) 



where (^k are random phases and the set of Vk's decreases 
as a given function of |k|, in agreement with experimental 
data |3- In principle, one should use for w(k) the disper- 
sion relation for electric drift waves (which are thought 
to be responsible for the observed turbulence) with a fre- 
quency broadening for each k in order to model the ex- 
perimentally observed spectrum ^(k, w). Unfortunately 
this would be prohibitive from a computational point of 
view, therefore one is led to simplify the model drasti- 
cally by choosing a;(k) = constant and the phases (/3k 
at random to reproduce a turbulent field (with the rea- 
sonable hope that the properties of the realization thus 
obtained are not significantly different from their aver- 
age). In addition we take for |Vk| a power law in |k| 
to reproduce the spatial spectral characteristics of the 
experimental 5(k), see Ref. 7]. Thus we consider the 
following explicit form for the electric potential 



V{x,v,t) = — 



N 



3in [^{nx + my) + (p^m - i^ot] 



(3) 

By rescaling space and time, we can always assume that 
L = 1 and = 1- In what follows, we choose N = 25. 
The spatial coordinates x and y play the role of the 
canonically conjugate variables. We extend the phase 
space {x,y) into {x,y,E,T) where the new dynamical 
variable r evolves as T{t) = t + T{0) and E is its canonical 
conjugate. We absorb the constant c/B of Eq. in the 
amplitude a, so that we can assume that a is small when 
B is large. The autonomous Hamiltonian of the model is 



Hix,y,E,T)=E + Vix,y,T) 



The equations of motion are 
dH dV dH 



dy dy 



y 



dx 



dV_ 

dx 



(4) 



1, (5) 



and E is given by taking H constant along the trajecto- 
ries. Thus, for small values of a, Hamiltonian Q is in the 
form H — Ha + eV, that is an integrable Hamiltonian Hq 
(with action-angle variables) plus a small perturbation 
eV. For simplicity we assume that the average of V over 
the angles is zero. Otherwise similar calculations follow- 
ing these lines can be done. The problem of control in 
Hamiltonian systems is to find a small perturbation term 



where {• ,•} is the Poisson bracket and the operator S 
is acting on H as SH — {S, H}. An expansion in power 
series in e of H'e gives 



H'e = Ho + e[{S,Ho} + V] 



f2 + {S,V} + -{S,{S,Ho}} 



The generating function S is chosen such that 
{S, Ho} + V = 0, 



(7) 



(8) 



provided that this equation has a solution. The control 
term /2 given by the cancellation of order terms [and 
using Eq. ©], 



f^ = --{S,V}, 



(9) 



satisfies the required condition that He is canonically con- 
jugate to Hq up to order terms. We notice that by 
adding higher order terms in e in the control term, one 
can build / such that He = H + e^f is integrable for 
sufficiently small e (see Ref. Q for more details and for 
a general formulation and theorem). 
In the case we consider, Hq = E and Eq. ^ becomes 



+ ^ = 



(10) 



and so S is one primitive in time of V. We choose the one 
with zero time average. For the model the generating 
function S is 



S{x,y,T) = — 

ZTT 



N 

E 

2 + m2<Jv2 



COS [27r(nx -I- my) + (p„ 



and the computation of /2 gives 



(11) 



/2(a;,y,r) 



E 



n2mi — nim2 



ni,mi,n2,m2 



m 



2)3/2( 



X sin[27r[(ni - n2)x + (mi - m2)y] + ipmmi - •^n^i-m^ ■ 

(12) 

We note that for the particular model (j3Jl, the partial 
control term /2 is independent of time. 

With the aid of numerical simulations (see Ref. 01 for 
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FIG. 1: Poincare surface of section of a trajectory obtained 
for Hamiltonian (|3J assuming a = 0.8. 
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FIG. 2: Poincare surface of section of a trajectory obtained for 
the same initial condition as in Fig.l and adding the control 
term (I12II to Hamiltonian ^ . 
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FIG. 3: Mean square displacement {r'^{t)) versus time t 
obtained for Hamiltonian Q with three different values of 
a = 0.7,0.8,0.9. 
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FIG. 4: Diffusion coefficient D versus a in log-log scales ob- 
tained for Hamiltonian @ (open squares) and Hamiltonian 
^ plus control term (full circles). 



more details on the numerics) , we check the effectiveness 
of the above control by comparing the diffusion prop- 
erties of particle trajectories obtained from Hamiltonian 
and from the same Hamiltonian with the control term 
(|12|l . Figures 1 and 2 show the Poincare surfaces of sec- 
tion of two trajectories issued from the same initial con- 
ditions computed without and with the control term re- 
spectively. Similar pictures are obtained for many other 
randomly chosen initial conditions. A clear evidence is 
found for a relevant reduction of the diffusion in presence 
of the control term 112|l . 

In order to study the diffusion properties of the sys- 
tem, we have considered a set of particles (of or- 
der 100) uniformly distributed at random in the domain 
< x,y < 1 at t = 0. We have computed the mean 
square displacement {r'^{t)) as a function of time 

1 ^ 

i— 1 

where Xi(i), z = 1, . . . , is the position of the i-th par- 
ticle at time t obtained by integrating Eq. |(SJ with initial 
condition Xj(0). Figure 3 shows {r'^{t)) for three different 
values of a. For the range of parameters we consider, the 



behavior of {r^{t)) is always found to be linear in time for 
t large enough. The corresponding diffusion coefficient is 
defined as 

..^limm. 

t^co t 

Figure 4 shows the values of I? as a function of a with 
and without control term. It clearly shows a significant 
decrease of the diffusion coefficient when the control term 
is added. As expected, the action of the control term 
gets weaker as a is increased towards the strongly chaotic 
phase. 

We check the robustness of the control by replacing 
/2 by 5 • /2 and varying the parameter 5 away from its 
reference value 5 = \. Figure 5 shows that increasing or 
decreasing 5 from 5 = \ result in a loss of efficiency of the 
control. The fact that a larger perturbation term (^ > 1) 
does not work better than the one with 5—1 means that 
the perturbation is "smart" and that it is not a "brute 
force" effect. 

Let us define the horizontal step size (resp. vertical 
step size) as the distance covered by the test particle 
between two successive sign reversals of the horizontal 
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FIG. 5: Diffusion coefRcient D versus S magnitude of the 
control term 1121 1 for a fixed value of a = 0.7 
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FIG. 6: PDF of the magnitude of the horizontal step size with 
and without the control term. 



(resp. vertical) component of the drift velocity. The ef- 
fect of the control is analyzed in terms of the Probabil- 
ity Distribution Function (PDF) of step sizes. Follow- 
ing test particle trajectories for a large number of initial 
conditions, with and without control, leads to the PDFs 
plotted in Fig. 6. A marked reduction of the PDF is ob- 
served at large step sizes with control relatively to the 
not controlled case. Conversely, an increase is found for 
the smaller step sizes. The control quenches the large 
steps, typically larger than 0.5. 

In order to measure the relative magnitude between 
Hamiltonian Q and /2, we have numerically computed 
their mean squared values: 



0.13a 



(14) 



which means that the control term can be considered a 
small perturbative term (when a < 1). 

So to conclude this work, we have provided an effective 
new strategy to control the chaotic diffusion in Hamil- 
tonian dynamics using small perturbations. Since the 
formula of the control term is explicit, we are able to 
compare the dynamics without and with control in a sim- 
plified model, describing anomalous electric transport in 
magnetized plasmas. Even though we use a rather simpli- 
fied model to describe chaotic transport of charged parti- 
cles in fusion plasmas, our result makes conceivable that 
to apply some smart perturbation can lead to a relevant 
reduction of the turbulent losses of energy and particles 
in tokamaks. 



[a] e-mail: ciraolo@arcetri.astro.it 

[b] e-mail: chandre@cpt.univ-mrs.fr 

[c] e-mail: lima@cpt.univ-mrs.fr 

[d] e-mail: pettini@arcetri.astro.it 

[e] e-mail: vittot@cpt.univ-mrs.fr 

[f] e-mail: charles.figarella@cea.fr 

[g] e-mail: philippe.ghendrih@cea.fr 

[1] G. Chen and X. Dong, Int. J. Bif. & Chaos, 3, 1363 (1993). 

[2] R. Lima and M. Pettini, Phys. Rev. A 41, 726 (1990); Int. 
J. Bif. & Chaos 8, 1675 (1998); L. Fronzoni, M. Giocondo, 
M. Pettini, Phys. Rev. A 43, 6483 (1991). 

[3] Interestlingly enough, the dynamics of test particles 
streaming along a field line in a fusion device are also gov- 
erned by Hamiltonian dynamics, the Hamiltonian being 
either the poloidal or toroidal magnetic flux. The princi- 



ple of control of stochastic transport would apply equally 
well to magnetic turbulence. 

[4] M. Pettini et ai, Phys. Rev. A 38, 344 (1988). 

[5] M. Pettini, "Low Dimensional Hamiltonian Models for 
Non-Collisional Diffusion of Charged Particles", in Non- 
Linear Dynamics (Bologna, 1988), G. Turchetti, Ed. 
World Scientific (Singapore, 1989) 287. 

[6] B. D. Scott, Phys. Plasmas, 10, 963 (2003) and the list of 
references [3] therein. 

[7] A.J. Wootton et al, Phys. Fluids B, 2, 2879 (1990). 

[8] M. Vittot: Perturbation Theory and Control in Classi- 
cal or Quantum Mechanics by an Inversion Formula , 
archived in jhttp : // arxiv . org/ pdf / math-ph/ 030305 1 1 



